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The coefficient of self-diffusion for a homogeneously cooling granular gas changes significantly if 
the impact- velocity dependence of the restitution coefficient e is taken into account. For the case of a 
constant e the particles spread logarithmically slow with time, whereas the velocity dependent coef- 
ficient yields a power law time-dependence. The impact of the difference in these time dependences 
on the properties of a freely cooling granular gas is discussed. 

PACS numbers: 81.05.Rm, 36.40. Sx, 51.20.+d, 66.30.Hs 
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I. INTRODUCTION 

The behavior of freely evolving granular gas has been 
intensively discussed recently, in particular the process of 
cluster formation due to inelastic collisions was of wide 
interest, e.g. [|]-|3|. On the basis of a continuum descrip- 
tion the effect of clustering in a force-free granular gas 
has been explained as an instability of the hydrodynamic 
equations [jl]-|| ■ F° r a deeper understanding of clustering 
phenomena it may be worth to consider the processes in 
a granular gas which precede clustering. To this end we 
investigate the effect of self-diffusion of particles in the 
regime of homogeneous cooling. 

The collisions of particles in a granular gas can be de- 
scribed by the coefficient of restitution e which relates the 
normal component of the relative velocity Vy = Vi — Vj 
before a collision to that after the collision as |t&e| = 

\j/\fij\ gives the direction of 



r 



e\vye\. The unit vector e 
fij = fi — fj at the instant of the collision. 

Up to our knowledge, all analytical calculations for the 
force-free case reported so far refer to systems of particles 
colliding with a constant restitution coefficient e. Experi- 
ments as well as theoretical studies show, however, that e 
depends on the normal component of the impact velocity 
\vijS\ §■ 

The problem of the restitution coefficient's dependence 
on the impact velocity has been addressed in where 
the generalization of the Hertz contact problem was de- 
veloped for the collision of viscoelastic particles. From 
this generalized Hertz equation one obtains the velocity- 
dependent restitution coefficient as an expansion: 
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where Y is the Young modulus, v is the Poisson ratio, 
and A depends on dissipative parameters of the parti- 
cle material (for details see ||). The effective mass and 
radius is defined as 



R eff ^ R 1 R 2 /(Ri + R 2 ) 
m e ff = m\mil (m,\ + 7112) 



(3) 
(4) 



with R1/2 and being the radii and masses of the 

colliding particles. The constants C\ = 1.15344 and 
C2 = 0.79826 were obtained analytically (6), confirmed 
then by numerical simulations and may be also written 
in a closed form as lf7|] 
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Eq.(|TJ) refers to the case of pure viscoelastic interaction, 
i.e. when the relative velocities \vijS\ are not too large 
(to avoid plastic deformation of the particles) and are 
not too small (to allow to neglect surface effects such 
as roughness, adhesion and van der Waals interactions). 
This implies that the initial temperature of the granular 
gas is not too large and the final temperature is not too 
small. The range of validity of Eq. (Q) depends on mate- 
rial and surface properties. In spite of such restrictions 
some important systems in nature do exist (e.g. plane- 
tary rings), where these conditions are satisfied j8j. Here 
we assume that the granular gas conditions allows for the 
application of Eq.(Q). 

For equilibrium 3D-system the time-dependence of the 
mean-square displacement reads 



{Ar{t)Y) =6Dt. 

cq 



(7) 



with 



where (• • •) denotes the equilibrium ensemble averaging. 
To calculate the mean-square displacement, one writes 



2 YVReff 
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(2) 
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and encounters then with the velocity autocorrelation 
function, (v(t')v(t")) , For systems at equilibrium this 
depends only on the time difference, |t' — t"\ and decays 
with a characteristic time t v . Therefore, at time t>r„ 
one arrives at the self-diffusion coefficient 
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(v(0)v(t) )„„ dt . 
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(9) 



For nonequilibrium systems such as granular materials 
the concept of the self-diffusion coefficient may be also 
applied, but with some care and with necessary general- 
ization. Obviously, this refers only to "liquid" or gaseous 
phases of the material || where the particles have a no- 
ticeable mobility. In the following we also restrict our- 
selves to homogeneous cooling and consider the stages of 
evolution preceding the cluster formation i.e. we 

assume that the granular material is (at least locally) 
homogeneous and isotropic. Hence, one can define the 
temperature T(t), which decreases with time due to the 
loss of energy according to inelastic collisions. For the 
impact-velocity dependent restitution coefficient ([!]) one 
has [||: 



T = To/(l+i/ro) 5/3 



(10) 



where To is the initial temperature and to is the charac- 
teristic time of the cooling process, which may be esti- 
mated as 
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(11) 



with a — 2R and n being the particle diameter and the 
particle number density, respectively. The mean-collision 
time 



- 1 {t)=ATi 1 ' 2 g 2 {a)na 2 T 1 ' 2 



(12) 



depends on time via the time-dependent temperature. 
Here g 2 (ct) is the contact value of the two-particle radial 
distribution function and the particles are of unit mass. 
Thus, the ratio of the two characteristic times reads 



M/ro ~ <5 11/6 [*/r c (0)] 



5/6 



(13) 



where r c (0) 1 is the collisional frequency at initial time, 
and S — Aa 2 / 5 T^ 10 is supposed to be small. Clearly, one 
can employ the concept of temperature if T c (t)/ro <C 1. 
Thus, Eq. (|l|) gives the estimate t < r c (0)5^ 11/5 for the 
upper time limit for which the use of the local tempera- 
ture is justified. 

An important property of homogeneously cooling gran- 
ular gas is that the velocity distribution is close to 
Maxwellian. Moreover, it persists with time, changing 
in accordance with the decreasing temperature ||. The 
small value of the fourth cumulant of the velocity distri- 
bution function for any value of the restitution coefficient 
reported in Jl^] also supports the Maxwellian distribu- 
tion p0 11 1. Therefore, we assume that the Maxwellian 



distribution and the molecular chaos hypothesis may be 
used with a good degree of accuracy. 

The evolution on the hydrodynamic time-scale may be 
described using the kinetic coefficients calculated on the 
short-time scale t ~ t c . For granular gases these trans- 
port coefficient will be now time-dependent. We calcu- 
late the self-diffusion coefficient within the uncorrelated 
binary collisions approximation and assume that the in- 
equality r c (t) <C r always holds true (see the above dis- 
cussion on this condition). Thus, on the time-scale t ~ r c 
the temperature may be considered as a constant. For 
t ^ r c , however, the self-diffusion coefficient occurs to be 
time-dependent and the generalization of Eq. (Q) reads 



^(Ar(<)) 2 ^ = 6 f D{t')dt' 



(14) 



where (• • •) denotes averaging over the nonequilibrium en- 
semble, which evolution is described by a time-dependent 
iV-particle distribution function p(t). Here we restrict 
ourselves to times when the hydrodynamic contribution 
to the self-diffusion coefficient is not large pUfl , so that 
( |l4| ) with D(t) calculated on the time scale t ~ t c gives 
an accurate description of the mean-square displacement. 

The aim of the present study is to analyze how the 
velocity dependence of the restitution coefficient influ- 
ences the diffusion in a gas of identical particles. The 
paper is organized as follows: In next Sec. II we obtain 
the time-dependence of the diffusion coefficient and the 
temperature of the granular gas in homogeneous cooling 
state. We also show that the mean-square displacement 
depends on time quite differently for the case of the con- 
stant restitution coefficient and for that determined by 
the impact velocity. In Sec. Ill we discuss our results for 
the mean-square displacement in a context of its possible 
impact on clustering. 



II. TIME DEPENDENCE OF THE DIFFUSION 
COEFFICIENT AND OF TEMPERATURE 



To describe the dynamics of the granular material we 
use the formalism of the pseudo-Liouville operator C |l4| ] 
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E^-^+E^ 



(15) 



i<j 



The first sum in ( |l5| ) refers to the free streaming of the 
particles (the ideal part) while the second sum refers to 
the particle interactions which are described by the bi- 
nary collision operators fl7| 



T, tj = a 2 jd 2 eQ 
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(16) 
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where Q(x) is the Heaviside function. The operator 6|- is 
defined as 



Hjf {n,rj,Vi,Vj •••) = / (n, fj,v-, v- ■ ■ ■) 



(17) 



where / is some function of dynamical variables. The 
after-collision velocities of the colliding particles, v[ and 



Vj are related to their pre-collisional values Vi, Vj via 



(18) 



The pseudo-Liouville operator gives the time derivative 
of any dynamical variable B (e.g. ]16[): 



B({fi,Vi},t) = iCB({r u Vi} ,t) 



(19) 



and, therefore, the time evolution of B reads (t > t') 

B ({*,*},*) = e^*-*'^ ({n,^},t') . (20) 
With (E0|) the time-correlation function reads 



(B(t')B(t)) = / tff>(f)B(* / )e a:( *~*' ) - B (* / ) , (21) 



where J" cfT denotes integration over all degrees of free- 
dom and p{i!) depends on temperature T, density n, etc., 
which change on a time-scale i>T c . In accordance with 
the molecular chaos assumption at t ~ r c the sequence 
of the successive collisions occurs without correlations. If 
the variable B does not depend on the positions of the 
particles, its time-correlation function reads ]lq] 

(B(t')B(t)) = (B 2 ) t , e -\ t - t '\' TB ^ (t > t') . (22) 

where (• • denotes the averaging with the distribution 
function taken at time t' . The relaxation time tb is 
inverse to the initial slope of the autocorrelation func- 
tion ]18[ |. It may be found from the time derivative of 
(B(t')B(t)) taken at t = t' . Eqs. <M) and @ then yield 



\t< 
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The relaxation time T B 1 (t / ) depending on time via the 
distribution function p(t'), changes on the time-scale 

t > T c . 

Let Bit) be the velocity of some particle, say V\it). 
Then with 3T(t) = (u 2 ) Eqs. (f|||) (with (|l||f)) read 



(^(t')-Ui(t))=3T(t / )e-l t - t 'l/ T ^ t ') 



-r l 7 1 (t') = (iV-l) 



(ui -Wi) t , 



(24) 



(25) 



To obtain ( |25| ) we take into account that £0^1 = 0, 
TjjUi = if i 7^ 1 and the identity of the particles. 



The calculation of r~ 1 (t') may be performed if we as- 
sume that the distribution function pit') is a product of 
the coordinate part, which corresponds to a uniform and 
isotropic system, and a velocity part which is a product 
of Maxwellian distribution functions 



m = exp[ ^ /2 I ( 5 )] . i = i. 



[2nT(t')} 3 / 2 

Integration over the coordinate part in (|25|) yields 



(26) 



(iV-1) p(t')S (fij-ae) dri ■ ■ ■ dr N =ng2(cr) , 

(27) 

where we use the definition of the configurational distri- 
bution functions Jll|, and where |l6| , |l8|| 



g2 (a) = ~(2-r ) )/il-r ] ) 3 



(28) 



gives the contact value of the configurational distribution 
function and i] — ^irna 3 . With 



(29) 



due to the collision rules and definition (16), one finally 
arrives at 



T v 00 = ^ng 2 (a)(7 2 j dv 12 <j){v 12 ) 



(30) 



d 2 e<d{-vi2 ■ e)\v\2 ■ e| (vi 2 ■ e) 2 (l + e)/ {v\ ■ v\) t , 



where 



(« 12 ) = (47rT)- 3 / 2 exp(- V 2 2 /4r) 



(31) 



is the Maxwellian distribution for the relative velocity of 
two particles. For e not depending on v\2 Eq. ( |30"| ) yields 

TvHt) = ^lna 2 g 2 (a)V*m = eJ ^Tz 1 it), (32) 

where Te(i) = § T c (t) is the Enskog relaxation time [fl6|| . 
For the granular gas it depends on time according to the 
same time-scale as the temperature. As shown in ( |32| ) the 
velocity correlation time for inelastic collisions exceeds 
that of elastic ones. This follows from partial suppression 
of the backscattering of particles due to inelastic losses 
in their normal relative motion. As a result the trajecto- 
ries of particles are more stretched as compared with the 
elastic case, and therefore the velocity correlation time is 
larger. 

As discussed above a constant restitution coefficient is 
not consistent with the nature of the inelastic collisions. 
Substituting (|l|) into ( pOf ) one finds the velocity correla- 
tion time for the gas of inelastically colliding spheres: 
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vio; 



dAa 2 / 5 (4T(i)) 1/10 - 



|r H± ) CfA 2 a^ (4T(t)) 1/5 T 



(33) 



where r(x) is the Gamma-function, te is given by 
( |32| ) pp[ | and we use (^|), which relates the coefficients 
C\ and C2. From (|3^) follows that the velocity autocor- 
relation function decays (as expected) on the short time 
scale, since t v is of the order of r c . 

Using the velocity correlation function one writes 



dt'3T(t') / dt"e-l*"- t 'l/ T »( i ') . (34) 



(Ar(f))' 



On the short-time scale t ~ r c , T(t') and Tt,(i') may 
be considered as constants. Integrating in (|34|) over t" 
and equating with ([l4|) yields for i ^> r c ~ r„ the time- 
dependent self-diffusion coefficient 



D(t)=T{t)r v (t). 



(35) 



Using the pseudo-Liouville operator one can also de- 
scribe the time-dependence of the temperature of the 
granular gas with the impact-velocity dependent resti- 
tution coefficient. From © it follows (see also §J|): 



Calculations similar to that for r v (t') yield 
T = -b 1 T 8 / 5 + b 2 T 17 / 10 T --- 



(36) 



(37) 



h = 4 • 2 1/5 7r 1/2 r 
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b 2 = ^2 2 /V/ 2 r (Ipj C 2 o 2 ng 2 {o)A 2 a^ 

Solving Eq. ([37]) and expanding the result in terms of the 
small parameter 



S = Aa 2/b T, 



2/5 T l/10 



(38) 



one arrives at 
T(t) (_ t 



To 



1- 
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t 



(39) 



with ai and a 2 being pure numbers pjl and with 
^-*-r.(0)-.| 4 -/-ftr(^) 



0.831928 • J • r c (0) _1 



(40) 



The leading term in this expansion corresponds to the 
dependence ( ]To| ) obtained previously using scaling ar- 
guments ||. From Eqs. (^,^,^3|) follows the time- 
dependence of the self-diffusion coefficient: 



+ a 3 (5 1+ — + a 4 <5 2 1+ — 
Do V W V T o J \ To, 



with pure numbers 03 and 04 |2l| and with 



n ° 1/2 2 1 \ rrl/2 



(41) 



(42) 



Correspondingly, the mean-square displacement reads 
asymptotically at To -C <: 

((Ar(<)) 2 )^i 1 /6 + a3(51ogt . (43) 

This dependence holds true for time 

r^O)^ 1 «t«r c (0),5- 11 / 5 , (44) 

where the first inequality follows from the condition 
To <C t, while the second one follows from the discussed 
condition T c (t) <C To. For the constant restitution coeffi- 
cient one obtains 



T(i)/To = [l + 7 ot/T c (0)]- 



(45) 



where 70 = (1 - e 2 )/6 @,||. Thus using Eqs.@ and 
(|35| ) one obtains for the mean-square displacement in this 
case 



((Ar(i)) 2 )~logi. 



(46) 



As it follows from Eqs.(|43|) and (|46|) the impact-velocity 
dependent restitution coefficient (pi) leads to a significant 
change of the long-time behavior of the mean-square dis- 
placement of particles in the laboratory-time. Compared 
to its logarithmically weak dependence for the constant 
restitution coefficient, the impact-velocity dependence of 
the restitution coefficient (|l|) gives rise to a considerably 
faster increase of this quantity with time, according to a 
power law. 



III. RESULTS AND DISCUSSION 

We studied the diffusion of particles in a homo- 
geneously cooling granular gas. With the assump- 
tion of molecular chaos we calculated the velocity 
time-correlation function and self-diffusion coefficient. 
For impact-velocity dependent restitution coefficient we 
found a relation which expresses the diffusion coefficient 
in terms of material constants of particles and character- 
istics of the granular gas, such as temperature, density, 
etc. 

In our calculations we used the restitution coefficient 
for viscoelastic collision of particles, which depends on 
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the normal component of the impact velocity |i%e| as 
a series £ |t%e|T with 7 = 1/5, 7 = 2/5,... [Eq.(g)]. 
Our approach can be extended to arbitrary exponents 
7, which let the expansion (Q) converge. Using these 
exponents one may possibly describe collisions with very 
large impact velocities (when plastic deformation occurs) 
or very small velocities (when the surface effects are im- 
portant), provided fragmentation/coagulation is ignored. 

For granular particles suffering viscoelastic collisions 
we found that the mean-square displacement grows with 
time as a power law ~ i 1 / 6 , i.e. much faster than the log- 
arithmic growth ~ logt, observed in granular gases with 
a constant restitution coefficient. It worth to note that 
qualitatively this power-law dependence (as well as the 
logarithmic one) simply follows from scaling arguments 
and the time-dependence of temperature. Indeed, the av- 
erage velocity scales as v ~ T 1 / 2 , and therefore as ~ t^ 1 
for the constant restitution coefficient and as ~ i~ 5 / 6 
for the impact- velocity dependent coefficient. The dif- 
fusion coefficient in granular gas scales as D ~ 1 2 /t c , 
where I ~ a^rC 1 is the mean- free path, which does 
not change with time (in the regime preceding cluster- 
ing), and r c ~ l/v is the mean-collision time. Thus, 
D ~ Iv ~ T 1 / 2 , and we obtain that the mean-square dis- 
placement, J D(t)dt, scales as ~ logt in the former case 
and as ~ i 1 / 6 in the latter case. 

What will be the impact of this apparently dramatic 
difference in the time dependence of ^(Ar(£)) 2 ^ on the 
properties of the granular gases? In the laboratory- 
time this corresponds to the enhanced spreading (and 
therefore mixing) of particles with the velocity-dependent 
restitution coefficient as compared with the case of con- 
stant e. Since the temperature decreases more slowly for 
the former case, as ~ i~ 5 / 3 , as compared with ~ i~ 2 for 
the latter one, retarded clustering may be expected. 

One might wish to define the average cumulative num- 
ber of collisions per particle Af(t) as a system inherent 
time-scale (t is the laboratory-time) and compare dynam- 
ics of the systems in their inherent-time scale. Af(t) is 
easily accessible in numerical simulations and is a conve- 
nient quantity to analyze evolution of granular gases (e.g. 
1 22,pi|). It may be found by integrating dJV = T c {t)~ 1 dt 
[ 13 1 . For constant restitution coefficient one obtains 
N(t) <~ logt, while for the impact-velocity dependent 
Af(t) ~ i 1 / 6 . Therefore, the mean-square displacement 
behaves identically in both cases as 



(Ar(AO) 2 ) - AT 



(47) 



If the number of collisions per particle Af(t) was the 
only quantity specifying the stage of the granular gas 
evolution, one would conjecture that the dynamical be- 
havior of a granular gas with a constant e and velocity- 
dependent e is identical provided jV-based time-scale is 
used. According to our understanding, however, this is 
not a adequate description of physical reality. Indeed, as 



it was shown in Ref. p2[ , the value of J\f c , correspond- 
ing to a crossover from the linear regime of evolution 
(which refers to the homogeneous cooling) to nonlinear 
regime (when clustering starts) may differ by orders of 
magnitude, depending on the restitution coefficient and 
on density of the granular gas. Therefore, to analyze 
the behavior of granular gas one can try an alternative 
inherent time-scale, T _1 = T(t)/To which is based on 
the gas temperature B . Given two systems of granular 
particles at the same density and the same initial tem- 
perature To, consisting of particles colliding with con- 
stant and velocity-dependent restitution coefficient, re- 
spectively, the time T allows to compare directly their 
evolution. A strong argument to use a temperature-based 
time has been given by Goldhirsch and Zanetti fl : They 
found that there are two main contributions to the rate 
of temperature decay. The first one refers to cooling due 
to inelastic collisions with a rate ~ T 3 / 2 (for constant e) 
and corresponds to the homogeneous cooling. The second 
one refers to viscous heating and behaves as ~ T 1 / 2 . Ini- 
tially the first is much larger than the second, but when 
temperature decays, the second contribution takes over 
and the shear wave adiabatically enslaves the tempera- 
ture field |[| . This corresponds to the nonlinear stage of 
evolution, when clustering starts. Thus, the temperature 
may indicate the stage of evolution of a granular gas. The 
recent numerical results of Ref. |E| strongly support our 
assumption: It was shown that while M c differs by more 
than a factor three for two different systems (Af c = 70 for 
a system with e = 0.9 and packing fraction cf> = 0.245, 
and N c = 23 for e = 0.6 and cj> = 0.05 Q) the values 
of T c , (defined, as T{JV C )/T Q ) are very close (T c w 0.0031 
for the first system and T c w 0.0027 for the second p2fl). 
These arguments show that one could consider T as a 
relevant time-scale to analyze the granular gas evolution. 

With T(M)/T ~ e ~ 2loAr for constant restitution co- 
efficient and T{M)/Tq ~ J\f~ 10 for the impact- velocity 

dependent one, we obtain ^(Ar) 2 ^) ~ logT for the case 
of constant e and 



(Ar)" 



7-1/10 



(48) 



for the velocity-dependent restitution coefficient. Thus, 
for the temperature-based inherent time-scale we again 
obtain a power-law dependence for the latter system 
and logarithmically- weak time-dependence for the former 
one. 

In conclusion, we found that that the impact-velocity 
dependence of the restitution coefficient e significantly 
influences the mean-square displacement of the particles 
in granular gas in laboratory-time. As compared with 
the logarithmically weak time-dependence found for a 
constant restitution coefficient, the impact-velocity de- 
pending coefficient (Q) yields a power-law (^) . It causes 
increasingly enhanced spreading of the particles through 
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the system. This together with the fact that the temper- 
ature decreases more slowly for the velocity-dependent 
e leads to the suggestion of retarded clustering in such 
systems. 

We also analyzed mean-square displacement of par- 
ticles using different inherent time-scales. We found, 
that while for the number-of-collision based time-scale 
both systems behave identically, for the temperature- 
based time-scale still power-law behavior is observed 
for the case of velocity-dependent coefficient and 
logarithmically-weak behavior for the case of constant 
e. 

We thank M. H. Ernst and I. Goldhirsch for valuable 
discussions. 
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